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Abstract. The elastic production of vector mesons at HERA is discussed from the 
theoretical point of view. We briefly review different models, their successes and short- 
comings. Main emphasis is put on recent issues in perturbative QCD calculations. Mod- 
els including the vector meson wave function are compared with an approach based on 
parton-hadron duality. We discuss several refinements of these models in some detail, 
including the important role of off-diagonal parton distributions. 

1 Introduction 

Why are we interested in elastic vector meson production? First of all the process 
7*p —> Vp provides us with well distinguishable experimental signals in a wide 
range of the j*p cm. energy W, the virtuality of the photon Q^, and the mass 
of the vector meson My ■ Quite some data are already available for V = p,4> and 
J/'f', and even for the heavy T first measurements were published recently.fj In 
the future the range in Q^ and W and the precision of the data will increase. 
This enables us to study vector meson production in detail in the very interesting 
regime where the transition from soft to hard QCD dynamics is expected (and 
already seen) to take place. In addition, there is hope to make use of the high 
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sensitivity of this process on the gluon distribution xg{x,Q ) in the proton to 
constrain this quantity at small values of x better than through other processes. 
In the following we first sketch the basic picture of elastic vector meson 
production. In Section 2 we briefly discuss different theoretical models which 
are not based on the two gluon exchange picture which is then introduced in 
Section 3. There, starting from the basic leading order result known for long 
time, we develop corrections which improve the leading order formula. In Section 
4 recent issues in pQCD calculations as off-diagonal parton distributions, the 
influence of the vector meson wave function and an alternative approach using 
parton-hadron duality are discussed. We mainly concentrate on diffractive p 
meson electroproduction, but the presented perturbative model is also successful 
in the case of J/'F and T. Section 4 contains our conclusions and outlook. 

1.1 The basic picture 

In Fig. n^ the basic picture for the process 7*p -^ Vp is shown: first the pho- 
ton with virtuality Q^ ~ —q^ fluctuates into a quark-antiquark pair. This qq 

^ For the discussion of experimental results on the production of light and heavy 
quarkonia see ||llj2,Hl and references therein. 
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fluctuation then interacts elastically with the proton p, where the zig-zag hue 
represents the (for the moment unspecified) elastic interaction with the proton. 
The 7*p centre-of-mass energy is denoted by W, 

W^ = {q+p)\ (1) 

whereas 

t^ip-pT (2) 

is the four-momentum transfer squared. (In the foUowing we wih mainly restrict 
ourselves to the case of small t.) The shaded blob at the right stands for the 
formation of the vector meson V, which, to leading order, has to form from the 
qq pair with invariant mass squared My. It is important to note that at high 



Y(q) 



t 
Fig. 1. Diagram for the elastic production of a vector meson V in 7*p collisions 



energy W corresponding to small values of x, 

Q^ + Ml 

the timescales involved in the problem are very different^: the typical lifetime 
of the 7* — > qq fluctuation as well as the time for the formation of the vector 
meson V is much longer than the duration of the interaction with the proton, 
i.e. Tj'^qq, Tqq^Y ^ T^. Therefore the basic amplitude factorizes, as sketched 
already in Fig. |l|, into the qq fluctuation, the interaction amplitude Aqq-^-p and 
the wave function of the vector meson V , 

AiYp^Vp) = i;^qq®Aqq+p®,P^ij, (4) 

and the process becomes calculable within various models]^ Formally it has been 
shown that for Q^ larger than all other mass scales in the process there is factor- 
ization into a hard scattering subprocess, non-perturbative (and off-diagonal, as 
will be discussed later) parton distributions and the meson wave function. [|5|] This 
strict proof of factorization holds for longitudinally polarized photons, whereas 
meson production through transversely polarized photons is shown to be sup- 
pressed by a power of Q. 

Let us now turn to the discussion of different models for the 7*15 interaction. 

^ This definition for x, which is often called ^ or xp, is common in diffractive physics 

and should not be confused with the ordinary Bjorken-a;, a;Bj = Q /{Q +W ). 
^ For a more detailed discussion of the ordering of the timescales see e.g. n . 
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2 Some non-perturbative models 

The following short section is far from being a review of this rich field, but is 
intended to give a hint at some non-perturbative models, which contrast the 
perturbative description of diffractive scattering, which is the main subject of 
this article. 

• We will not cover approaches based on vector meson dominance (see e.g. ^). 

• For Regge-phenomenology-based models of (one or two) Pomeron exchange we 
refer the reader to Q]. 

• The model of the stochastic QCD vacuum 

Dosch, Gusset, Kulzinger and Pirner [gj have developped a model of the interac- 
tion with the proton, which is similar to the semi-classical model of Buchmiiller 
discussed in [gj in the context of inclusive diffractive DIS. This model, originally 
used for hadron-hadron scattering, leads to linear confinement and predicts a 
dependence of the high-energy scattering on the hadron size. It gives a unified 
description of low energy and soft high-energy scattering phenomena. Dosch et 
al. approximate the slowly varying infrared modes of the gluon field of the pro- 
ton by a stochastic process. Via a path integral method they average over all 
possible field configurations. For the splitting of the photon into the qq pair 
and for the description of the vector meson they use light cone wavef unctions. 
Within their model they are able to calculate the Q^ dependence of the cross 
section, as well as the dependence on the momentum transfer t, da/dt, and the 
ratio of the longitudinal to the transverse cross section, L/T, where longitudinal 
and transverse refer to the polarization of the photon. Their results are in fair 
agreement with experimental data. There is no prediction for the W dependence 
of the cross section. 

• Rueter has extended the model of Dosch et al. to also describe the W de- 
pendence of the cross section. [|o| He achieves this by using a phenomenological 
model based on the exchange of one soft and one hard Pomeron, each being a 
simple pole in the complex angular momentum plane, similar to the Donnachie- 
Landshoff model W . For the very hard components of the photon fluctuations 
he treats the interaction perturbatively and achieves a good description of the 
experimentally observed transition from the soft to the hard regime. 



3 The two gluon exchange model 

To leading order in QCD the zig-zag line in Fig. |l], which stands for the elastic 
scattering via the exchange of a colourless object with the quantum numbers of 
the vacuum, can be described by two gluons. If the scale governing the (trans- 
verse) size of the photon fluctuation is large compared to the typical scale of 
non-perturbative strong interactions, i.e. if 

Q'»^QCD or M^:^AlcB, (5) 
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then the coupling of the two gluons to the qq fluctuation can be treated reliably 
within perturbative QCD (pQCD). Another kinematic regime, where pQCD is 
applicable, is high-i diffraction. There the hard scale which is needed to ensure 
the validity of the perturbative treatment is given by the large value of the 
momentum transfer i, and one expects high-i diffractin to be a good place to 
search for the perturbative Pomeron. n% 

It has been shown some time ago that due to the factorization property of the 
process the coupling of the two gluons to the proton can, in the leading logarith- 
mic approximation, be identified with the ordinary (diagonal) gluon distribution 
in the proton. ||l^ , |l3|jl4|] We will come back to this point later when discussing 
the importance of off-diagonal gluon distributions. 

3.1 The basic formula 

The basic leading order formula for diffractive vector meson production is given 

by lira 



^Cr rVA43^3 „, fn'^\2 



V 



d^(^*^-^^)..o = ^l8^^^i^'^(^'^)J ^ + ^)' (^) 



where a is the electromagnetic coupling and the gluon distribution is sampled 
at the effective scale 

q' = (q2 + a4)/4. (7) 

In Eq. (g) the non-relativistic approximation for the vector meson wave function 
is used and the coupling of the vector meson to the photon is encoded in the 
electronic width Pj^. Note that Eq. (^ is valid for t = 0. In the approach 
discussed in the following there is no prediction for the t dependence of the 
cross section, which is assumed to be of the exponential form exp(— 6|i|) with 
an experimentally measured slope-parameter 6, which may depend on the vector 
meson V and on Q^. On the other hand, Eq. (||) makes predictions for both the 
Q^ and the W dependence of the cross section for longitudinally and transversely 
polarized photons for all sorts of vector mesons, as long as either Q^ or My is 
large enough to act as the hard scale. It is obvious that the W dependence comes 
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entirely from the gluon distribution xg{x, Q ), which enters quadratically in the 
cross section. 



3.2 Improvements beyond the leading order 

In the following we will discuss several improvements of the leading order for- 
mula.^ 

• Eq. (@) contains only the leading imaginary part of the positive-signature 
amplitude 

A(xt{x-^ + {-x)-^) . (8) 



For more detailed discussions see e.g. [[L5lJlq] or the recent review ll 
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The real part of the amplitude can be restored using dispersion relations: 

Re^ = tan(7rA/2) IniA , (9) 

where A is given by the logarithmic derivative 

_ dlogA 
^ - 91og(l/:r) ■ ^^"^ 

For the case of p production, the contributions from the real part are roughly 
15%. For J/tf' production in the HERA regime they amount to approximateley 
20% and are even bigger for T production [|18Hl9[ , where larger values of x are 
probed. 

• In Fig. S one of four leading order diagramaj for the two gluon exchange 
model is shown with some kinematic variables which will be used below. In the 
general case the two giuons gi, g2 have different x, x' and transverse momenta 
It, ^'t- The leading logarithmic approximation of the i"^ loop integral (indicated 
by the circle in Fig. 0) leads to the identification with the integrated gluon dis- 
tribution xg{x,Q ) at the effective scale Q defined in Eq. (^. Beyond leading 
logarithmic accuracy one has to perform the £^ integral over the unintegrated 
gluon distribution f{x,i'^). This can lead to numerical results which are, de- 
pending on the kinematical regime, twice as big as the result from Eq. (@).||lE 




P ■ I ■ - P 
I 

Fig. 2. One of four leading order diagrams for the two gluon exchange model for 
diffractive vector meson production 



Although here we are considering elastic production at small momentum 
transfer t, the timelike vector meson with mass My has to be produced from the 
spacelike (or real) photon with virtuality Q^. This means, that even if there 
is no transverse momentum transfer, £t = £'tj there has to be a difference 
x — x' = {My + Q^) / (W^ + Q^) in the longitudinal momentum of the two 
giuons gi and g2- Therefore the identification with the ordinary diagonal gluon 

^ There are three similar diagrams: one where both giuons couple to the antiquark 
and two where one gluon is attached to the quark, whereas the other couples to the 
antiquark. 
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distribution x g{x,Q ) is only a good approximation for very small values of x 
and i, and in the general case the process 7*p -^ Vp depends on ofF-diagonal par- 
ton distributions. pQ] Their importance for diffractive vector meson production 
will be discussed in the following. 



4 Recent issues in pQCD calculations 

4.1 Off-diagonal parton distributions 

Off-diagonal (also called "skewed" or non-forward) parton distributions^ are 
much studied recently]^ In the case of small t scattering the skewedness comes 
from the difference between x and x' of the two gluons gi and g2, and the 
cross section can be shown to be proportional to the square of a skewed gluon 
distribution, 



c'g(a 



2 



(11) 



Here x = (M^^- + Q^) / {W^ + Q^) , x' = (A/|j- - M^) / {W^ + Q^) « x, 
and Mgg is the mass squared of the intermediate qq pair. (Taking the leading 
imaginary part of the amplitude corresponds to cutting the amplitude as indi- 
cated by the dashed line in Fig. g and putting both q and q on-shell, which in 
turn fixes x. x' has to accomodate the difference between Mqq and My and it not 
fixed due to the integration over all possible quark (and antiquark) momenta. 
At leading logarithmic order x' <S^ x and we can put x' ~ 0.) 

For arbitrary kinematics skewed parton distributions are not connected with 
the diagonal ones and are unknown non-perturbative objects. However, in the 
case of small x, they are determined completely by the diagonal ones.|2|,g2l The 
ratio of skewed to diagonal gluon distribution is given by 

_ x'g{x,x') 2^^+3r(A+|) 

xg{x) - ^ r{X + A)- ^'^^ 

Here F is the usual Gamma-function and the effective power A can be obtained 
from the logarithmic derivative of the amplitude A for the 7*p -^ Vp cross 
section, 

alog(l/x) 



As will be shown below, the magnitude of the resulting correction factor for the 

2 



total cross section, R^, can be sizeable, especially for large Q^ or M^ 



These off-diagonal parton distributions are not parton densities in the ordinary prob- 
abilistic sense but matrix elements of parton-fields between different initial and final 
proton states. 
^ See teol and references therein. 
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4.2 The vector meson wave function 

Another important issue is the treatment of the vector meson wave function. As 
sketched in Fig. |l| and Eq. (Q) it enters the amphtude via a convolution with 
the scattered qq fluctuation. In Eq. (|g) the non-relativistic approximation was 
adopted. This means, that quark and antiquark equally share the longitudinal 
momentum of the photon, i.e. z = I — z = 1/2, and that there is no inter- 
nal (transverse) momentum kx in the qq bound state. Therefore, in this naive 
approximation, 

i^^,iz,kT) = 5('\kT) 6 (z- 1/2), (14) 

and My = 2mq. While this simplification may be suitable for heavy mesons like 
the T, it is clear that the non-relativistic approximation has to break down for 
light quarks. Various groups have worked on improving this approximation by 
including the Fermi motion of the quarks in the meson by using a nontrivial 
wave function. |ll3| , |l6yiq , |2^ Different models for the meson wave function were 
used which lead to quite different correction factors: whereas in Gaussian models 
there is no strong suppression [^, the large fcy tail typical for wave functions 
from non-relativistic potential models seems to lead to large corrections fig] . 
On the other hand, considering that within these potential models a big part of 
the C'(w^) corrections comes from a regime, where kx is bigger than the quark 
mass itself, these large corrections may well be an artefact of the non-relativistic 
approximation. 

Another related problem is the question, which mass for the quarks should 
be used in the perturbative formulae. Note that Eq. ( |q) i s written in terms of 
the vector meson mass Aly- However, as discussed in |15| , the full expressions 
used to include higher order (relativistic) corrections contain the quark mass uiq 
instead of My- As the ratio My / {2mq) enters with a high power, this difference 
is not negligible and should be taken into account in the calculation of the 0{v^) 
corrections applied to Eq. (|^). 

In addition, it is well known that there are other relativistic corrections, 
which in principle have to be taken into account in a consistent way. As pointed 
out by Hoodbhoy |2^ , gauge invariance is only preserved if higher Fock states 
{qqg, qqgg, ■ ■ ■) are included in the wave function. In doing so he arrives at the 
conclusion, that the relativistic corrections to the quark propagators plus the 
corrections from the higher Fock states amount to only a few percent for J/'P 
production, in agreement with |l5|| . 

After all large relativistic corrections can probably be excluded, but, as dif- 
ferent approaches lead to quite different results there remains a considerable 
uncertainty and the issue a hot topic. 

4.3 An alternative approach based on parton-hadron duality 

In this section we will discuss an alternative approach, which avoids the meson 
wave function and leads to results which are in surprisingly good agreement with 
available data. It was proposed in p3 for p meson electroproduction, where the 
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hard scale is provided by Q^, not by Mpn Due to the tiny u and d quark masses, 
in the case of the p non-relativistic approximations cannot be justified, and the 
wave function is not very well known. Now the crucial problem was, that all 
naive predictions for the ratio of the longitudinal to the transverse cross section, 
which are based on the perturbative formula m), lead to 

at/^T - Q'A^p ■ (15) 

This is much too steep and incompatible with experimental data (see below). 
The inclusion of effects from a light cone wave function for the p does not change 
the picture considerably.n These observations indicate that the main effects are 
not coming from the p wave function and lead to the proposal of a different model 
in IgJ : there the cross section for p production is predicted via perturbative uu 
and dd quark pair electroproduction together with the principle of parton-hadron 
duahty (PHD) |g5|. PHD means that the integral of the parton (qq) production 
cross section over a mass interval AM is approximately equal to the sum over 
all (corresponding) possible hadron production cross sections in the same mass 
interval. In the region M^, w M^ the production of more complicated partonic 
configurations (like qq + g, qq + 2^, qq + qq, etc.) is heavily suppressed. On the 
hadronic side the p resonance (plus the small admixture of the tu) with its decay 
into two (three) pions completely saturates the cross section. Therefore we can 
well approximate the p production cross section 

7*p -^ pp -^ TTTrp 

by 

f * ^ no V^ r^' da (7> ^ iqq)p) 
aij p^pp)~ 0.9 2^ / -—5 (16) 



M2 



where Ma and Mb have to be chosen to embrace the p resonance appropriately, 
i.e. M^ — Ma ~ 1 GeV'^. The factor 0.9 on the right side of Eq. (|l6|) corrects for 
the contributions from u production. 

The perturbative formulae for the qq production cross section are derived 
from the amplitudes depicted in Fig. and can be written in terms of the con- 
ventional spin rotation matrices d'l (6) (see p3 for details): 



dVi 



dM^dt 
(Par 



dM^dt 



j dcose \d\^{e)\^\lL\\ (17) 

A'^'^ela M^ I f^ /, , ,2,1 ,2\ 9 



dcos0 {\d^Ae)r + \d^_je)r] \It 



3 (Q2 + M2)2 8 



t=o 3 (Q2 + M2)2 4 



Experimentally both the t dependence da /At ~ exp( — 6t) with & ~ 5 — 6 GeV~ for 
Q^ > 10 GeV^ and the W behaviour of the cross section a oc W'^'^ indicate that p 
meson electroproduction is not a soft, but mainly a hard process. 
^ One might argue that ax receives large contributions from the small kx region, which 
is non-perturbative. But those contributions would cause the transverse cross section 
to fall off even faster with increasing Q^ and therefore worsen the problem |24|. 
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where Cq is the electric charge of the quark g, a the electromagnetic coupling and 
6 the polar angle of the quark q in the qq rest frame with respect to the proton 
direction (fc-r = M/2sin0). Il_t are integrals over the gluon ^^ and given by 

Il{K') =K' J ^ aM) /(.,4) {-^2-^)' (18) 

with / being the unintegrated gluon distribution and 



K} = J{K-^+il^Y - 4^24, K^ = k^{Q'' + M^)/M\ 



In Eqs. (O) the different rotation matrices appropriately reflect the different 
spin states of the qq produced from longitudinal and transverse photons, and the 
integrals Il,t contain the scattering off the proton via the two gluon exchange .P^ 
In order to pick up only those uu, dd configurations which correspond to the 
quantum numbers of the p, one has to project out the J^ = 1 states. This 
can be easily done on amplitude level with the same rotation matrices d'l (6), 
see pj]. (Even higher spin states like the p{5~) can be projected out using 
the corresponding d function pa.) It is important to note that through the 
projection on amplitude level the longitudinal and transverse cross sections aL,T 
are less infrared sensitive than Eqs. (U^), and therefore ctt becomes calculable 
without a large uncertainty from the treatment of the (non-perturbative) infrared 
region. 

For the complete numerical predictions one also has to include the contribu- 
tions from the real part of the amplitudes and the skewed gluon distribution as 
discussed above. Both effects are taken into account on amplitude level. Eqs. ( [l7|) 
give the cross section differential in t for i = 0. To arrive at the t integrated total 
cross section one assumes the exponential form exp(— 6|i|). The slope b can be 
taken from experiment or theoretical models and depends in general on M'^, W 
and Q^ . For more details we refer the reader to |2j . 

To go beyond the leading order prediction in a completely consistent way 
would require in addition the full set of next-to-leading order ghionic correc- 
tions to the (qq)-2g vertex. These corrections are not known yetFj, but can be 
estimated by a /C factor [p|j24|. Similar to the Drell-Yan process, there are tt^ 
enhanced terms, which come from the in terms in the double logarithmic Su- 
dakov form factor. Resummation of those leading corrections results in the /C 
factor /C = exp {irCFCts) which leads to a considerable enhancement of the cross 
section. 



^^ Here we assume s channel helicity conservation (SCHC), i.e. the produced p has 
the helicity of the virtual photon. However, there are small violations of SCHC 
(e.g. 7y —> pl) which can be successfully described in a framework similar to the 
one discussed here. For a discussion of recent measurements of the 15 spin density 
matrix elements of p production compared to different theoretical predictions see_ 

^^ A first step towards the calculation of the full NLO corrections is provided by 
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35 40 45 
Q^ [GeV^] 



Fig. 3. a (7*p -^ pp) predicted within the PHD model as described in the text com- 
pared to recent HI data |P9|. Continuous line: using the skewed gluon distribution, 
dashed line: without skewing 



In Fig. the complete numerical prediction for ^*p — > pp using the PHD 
mode]0 is shown as a function of Q^ together with recent HI data.|22] The 
continuous line includes all the effects discussed above, whereas the dashed line 
does not include the skewed gluon. The importance of the off-diagonal gluon 
for the Q^ behaviour of the cross section is obvious and the effect seems to be 
required to describe the data. Of course the model prediction is not free from 
uncertainties like the choice of the mass interval M"^ — A/^ in Eqs. (|6|) or the 
scale of as in the /C factor. These (and other) uncertainties are discussed in detail 
in ||27| , but they affect mainly the normalization of the cross section and do not 
spoil the good agreement with the experimental data. In Fig. ^ the prediction of 
the PHD model for the ratio L/T is shown as a continuous line. It agrees fairly 
well with the data points, which show a very modest rise with Q^ in contrast to 
the naive prediction Eq. ( [l5| ) (dashed line). Thus, in the PHD picture, it is not 
the p wave function, but the dynamics of the qq pair creation from longitudinal 
and transverse photons together with the off-diagonal two gluon interaction and 
the projection onto the 1~ state, that determines the Q^ dependence and the 
ratio L/T. 

It is important to note that the PHD model also works in the case of massive 
quarks and heavy mesons. Starting from formulae for diffractive heavy quark 
production Q and modifying the projection formalism appropriately, elastic 
T photoproduction was recently predicted using PHD in agreement with first 
measurements. [n8| The same formalism can also be applied to diffractive J/!/' 
production, [g^ Again, as shown in Fig. ||, there is a surprisingly good agreement 
between the predicted cross section as a function of Q^ and the experimental 
data isl. 



For the numerical analysis the MRST99 gluon po| was used and the scale of Qs in 
the K. factor was chosen as 2K^ . For more details see u'm. 
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Fig. 4. L/T predicted within the PHD model {continuous line) as described in the text 
compared to experimental data P9| , pl| , p2[ and the naive prediction (dashed line) from 
Eq. di) 



first point is for Q" = 




Q^ [GeV^] 

Fig. 5. Cross section for diffractive J/'I' production as predicted in the PHD model 
127]] compared to recent HI data |33| 



5 Summary 

Elastic vector meson production is a rich field, both from experimental and the- 
oretical point of view. Different theoretical models describe the data, and more 
precise data in an increased kinematical range will be needed to clarify the sit- 
uation. We have briefly discussed some non-perturbative models, but mainly 
concentrated on perturbative approaches. We have shown that with recent im- 
provements pQCD-based approaches work very well and are in agreement with 
data. The fairly large impact of skewed parton distributions on the predictions 
within these models is supported by the data. There is good hope that in future 
we will be able to discriminate between the different models and to understand 
elastic vector meson production in more detail. By combining different observ- 
ables from different processes elastic vector meson production with its high sen- 
sitivity to the gluon at small x will finally help to constrain the gluon much 
better. For this much effort will be needed also from the theoretical side in order 
to increase the precision of the calculations. 
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